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Abstract 

A new gauge (m-gauge) condition is proposed by means of a gen- 
eralization of the Maximal Abelian gauge (MAG). The new gauge 
admits a space time dependent embedding of the residual U(l) into 
the SU(2) gauge group. This embedding is characterized by a color 
vector rh{x). It turns out that this vector only depends of gauge in- 
variant parts of the link configurations. Our numerical results show 
color ferromagnetic correlations of the rh{x) field in space-time. The 
correlation length scales towards the continuum limit. For compar- 
ison with the MAG, we introduce a class of gauges which smoothly 
interpolates between the MAG and the m-gauge. For a wide range of 
the gauge parameter, the vacuum decomposes into regions of aligned 
vectors rh. The "neutral particle problem" of MAG is addressed in 
the context of the new gauge class. 
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1 Introduction 



Revealing the mechanism which operates in Yang-Mills theories for confining 
quarks to color singlet particles is one of the most important tasks of modern 
quantum field theory. A knowledge of this mechanism will help to understand 
nuclear forces from first principles and will, e.g., have a strong impetus on 
the present understanding of hadron physics. 

In the recent twenty years of active research, many proposals have been 
launched to explain quark confinement (see [1-7] for an incomplete list). 
Nowadays large scale computer simulations assign top priority to the so- 
called color superconductor mechanism 0. This mechanism applies in the 
Abelian gauges J7J, which allow for a residual U(l) gauge degree of freedom. 
Projection is performed to reduce the full SU(2) to compact U(l) gauge the- 
ory. After this projection, monopoles which carry quantized color-magnetic 
charges with respect to the residual U(l) group naturally appear as degrees 
of freedom ||. The color-superconductor mechanism operates as follows: a 
condensation of these monopoles implies a (dual) Meissner effect. Color- 
electric flux is squeezed into flux tubes implying that the potential between 
two static color sources is linearly rising at large distances. Modern com- 
puter facilities provide the testing grounds: evidence has been accumulated 
that the dual superconductor picture captures parts of the roots of quark 
confinement ||, |K|. 

Unfortunately, the color-superconductor picture suffers from the following 
drawback: color states which are neutral with respect to the residual U(l) 
gauge are insensitive to the condensate of the color magnetic monopoles and 



therefore do not acquire a confining potential [11]. In this case, one would 
expect additional "light" states in the particle spectrum besides hadrons 
and glueballs. This contradiction to the experiment requires a refinement 



of the dual superconductor picture. For this purpose, it was argued [TT 
that all color-magnetic monopoles which are defined by different Abelian 
projection schemes condense while only those monopoles which correspond 
to the gauge fixing at hand are manifest. The idea of the condensation of 
"hidden" monopole degrees of freedom might conceptually solve the "neutral 
particle problem" , but conceals the non- Abelian nature of the superconductor 
mechanism. 

In this paper, we throw a new glance onto the non-Abelian Meissner effect 
by refraining from a residual U(l) gauge group which is uniquely embedded 
into the SU(2) gauge group all over space time. Instead of, space-time is 
decomposed into regions in each of which the embedding of the residual 
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U(l) group into the SU(2) gauge group is chosen to yield the minimal error 
by a subsequent projection. Particles which are "neutral" with respect to 
the U(l) subgroup in one particular region carry charge in another region. 
If the average over all configurations is performed during the Monte-Carlo 
sampling, all particles are confined on distances which exceed the intrinsic 
length scales of the regions. 

For putting this idea on solid grounds, we propose a new type of gauge 
(referred to as m-gauge) which appears as a generalization of the Maximal 
Abelian gauge (MAG). In the m-gauge, the orientation of the U(l) subgroup 
in SU(2) is specified by a unit-color vector fh which depends on space-time. In 
a first step, we show by means of numerical calculations that the m-projected 
theory still bears quark confinement at full strength. Secondly, for putting 
the m-gauge into a proper context, we investigate, by virtue of a gauge fixing 
parameter, a gauge fixing which smoothly interpolates between the MAG and 
the m-gauge. For a wide range of the gauge fixing parameters, the vacuum 
decomposes into regions of (approximately) aligned color vectors fh. 

The paper is organized as follows: in the next section, we will introduce the 
novel type of gauge fixing and the corresponding projection. The numeri- 
cal errors induced by projection are studied for the case of MAG and the 
case of m-gauge, respectively. In section [| we reveal the color ferromag- 
netic correlations of the vector fh present in the m-gauge. The interpolating 
gauges are introduced and the vacuum structure in these gauges is discussed. 
Conclusions are left to the final section. 

2 M-gauge 

2.1 The new gauge fixing and projection 

A particular configuration of the SU(2) lattice gauge theory is represented 
by a set of link matrices U^x) G SU(2), which are transformed by a gauge 
transformation Q(x) G SU(2) into 

-> U%{x) = Q(x) U^x) Sl\x + ft) . (1) 

Our new proposal for the gauge fixing condition is given by 




maximum , (2) 



m(x) = m a (x) r a , fh T (x)fh(x) = 1 , 
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where r a are the SU(2) Pauli matrices. For a given configuration U^(x) we 
allow for a variation of the gauge matrices Q(x) and of the auxiliary unit 
vector m a (x) for maximizing the functional Sfi x . Note that a reflection of 
the vector, i.e., m(x) — > —m(x), does not change the gauge fixing functional. 
Identifying the points rh and — rh of the sphere S2 defines a projective space 
RP2 which carries the gauge fixing information. Furthermore, is invariant 
under a multiplication of the gauge matrix Q, under consideration, with a 
center element of the SU(2) gauge group, i.e., Q — > (-l)Q. This implies 
that the theory after gauge fixing possesses a residual Z 2 gauge invariance, 
at least. It turns out that a further invariance is unlikely to exist for generic 
link configurations U^x) (see discussion in subsection |2~2j ). 

The concept of projection is to reduce the number of degrees of freedom while 
preserving the confining properties of SU(2) gauge theory. It might be easier 
in the reduced theory to reveal the mechanism of confinement than resorting 
to the full SU(2) gauge theory. In the present case, we define the projected 
links Un(x) by 



N 



U^{x) + m(x)U^{x)m(x) 



(3) 



where the normalization N is obtained by demanding U^x) UU 



x 



1. 



and the projection 
1 Vfi,x (4) 



It is convenient for an illustration of the gauge fixing 
(IH) to decompose the link variable as 

== a^\x) + ia,(x)r, (aj?>(s)) 2 + Ufa) 

In this case, the gauge fixing condition Sex (@) becomes 

s„ = E {W^)) 2 

fi,{x} 

Representing the vector a^x) by components parallel to rh(x) and perpen- 
dicular to m(x), i.e., a^(x) = (aj[, a 1 ^-, aj l ± ) T , the condition @ is equivalent 
to 

,2 / „,, ,\2 



a M (x) + 2 (maj 



maximum . 



(5) 
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fix 



E 

li, {x} 



1 



maximum . 



(6) 



This equation tells us that the gauge fixing 
a M (x) perpendicular to the vector m(x). 



minimizes the link components 



Inserting (f|) in 



one finds for the projected link variables 
= 2N 



(7) 



2N 



(x) + i aMx) fh(x) t 
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Projecting link configurations is a two-step process: firstly, one exploits the 
gauge degrees of freedom to minimize the link components a* (a;), aV~(x), 
perpendicular to m(x), and, secondly, these components ah L (x), aV~{x) are 
dropped for obtaining the projected link variable U^x). For quantifying the 
error of this projection, we introduce 

^ := it (S%T[U]) v , (9) 

-'Vlink 

where A cinic is the number of lattice links and S™ X [U] is the maximum value of 
the gauge fixing functional Sfi x (H) for a given link configuration U^(x). The 
brackets in (Q) denote the Monte-Carlo average over all link configurations. 
An inspection of the equations @ and @ tells us that projection yields the 
exact result if u possesses the largest value possible, i.e., u = 1. The error 
increases if u decreases. 

Note that the space-time dependence of m(x) relative to the link variables 
U^(x) in the functional Sr x (§) is obtained by the demand for minimizing 
the error induced by projection. This demand dictates that m(x) cannot 
be identified with a so-called Higgs field which figures in general Abelian 
gauges fll^| . For an illustration of this fact, let m(x), Q(x) be the configura- 
tions which maximize Sfi x [i7] (@) for a given link configuration [/^(x), and let 
assume that the fields m(x), Q(x) are uniquely defined (this is the generic 
case; see subsection [2.2| ). We introduce U^(x) as the link variables which 
are obtained from U^x) by the gauge transformation V(x), and repeat the 
gauge fixing procedure with the U^(x) as basis. If m v (x), Q v (x) denote 
the configurations which correspond to the maximum of S^U ], one finds 
fl v (x) = Q(x) V(x) and m v (x) = m(x). The later relation tells us that m(x) 
only depends on the gauge invariant parts of the link variables U^x), and 
therefore encodes physical information. By contrast, an auxiliary Higgs fields 
transforms homogeneously under the gauge transformation V{x). 

Let us compare our new gauge, defined by (0), with the Maximal Abelian 
gauge (MAG). The latter gauge can be obtained from the gauge condition (|2|) 
if one does not allow a variation of m(x) with space time. For constant vectors 
m, one might choose rh to point in three direction in color space without a loss 
of generality. For an SU(2) gauge theory in four dimensions, one associates 
four link variables with each space time point, and therefore counts 12 degrees 
of freedom at each lattice site (9 physical and 3 gauge degrees of freedom). 
The MAG projection effectively reduces the SU(2) gauge theory to an U(l) 
one. After projection, the number of degrees of freedom per site is therefore 
4. In the new gauge, presented here, naive counting yields four Abelian links 
and the unit vector m(x), i.e., 6 degrees of freedom at each site. 
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We finally comment on the Gribov problem to round out this subsection. 
Note that the following remarks also apply to the class of general Abelian 
gauges, and that, in particular, the practical problem in implementing the 
gauge is not a specialty of the new gauge proposed in the present letter. 
Although the gauge fixing starting from the condition (|2|) is conceptually free 
of Gribov ambiguities if one seeks out the global maximum of the functional 
Sfix (@) , one recovers the Gribov problem in practice when the algorithm fails 
to detect the global maximum. In the context of variational gauges, several 
strategies have been proposed for evading this problem. One possibility is 
introducing a Laplacian version of the gauge fixing condition for adapting the 
problem to the numerical capabilities. Results are available in the literature 
for the case of Landau gauge JI3|| , for the case of MAG [H| and for the 



case of the center gauge flSfl . Another possibility is to introduce quantum 



gauge fixing |TB| for putting the gauge fixing which is implemented by the 
algorithm in the proper context. Here, we will not perform a detailed study 
of the "practical" Gribov problem in the context of the new gauge ([|). In the 
present first investigation, we will only check whether the numerical results 
(see next section) are stable against random gauge transformations on the 
link variables before invoking the gauge fixing algorithm. 



2.2 Quality of projection 

Our numerical simulations were performed using the Wilson action and a 
lattice with 12 4 space-time points. For /5-values in the scaling window, i.e., 
(3 G [2.1,2.5], 200 heat-bath steps Jl7| were performed for initialization. 
When gauge fixing is requested by the application of interest, we used a 
standard iterative procedure with over-relaxation for finding the maximum 
value of the gauge fixing functional S^ x (0). Once we have obtained con- 
figurations m(x), Q(x) which maximize the functional Sfi x , we distort these 
configurations and re-handle the gauge fixing. Repeating this procedure (for 
a given link configuration U^(x)) several times, we find unique fields m(x), 
Q(x) at the maximum of Sfi x . This provides numerical evidence that the 
(local) maximum of S& x is stable against small gauge transformations, and 
that flat directions in the configuration space of m(x), Q(x) do not exist for 
generic link configurations. A thorough study of the Faddeev-Popov deter- 
minant in the case of m- gauge is requested for a rigorous proof of this fact. 
This is left to future work. 

In a first investigation, we calculated the Creutz ratios \kk with help of the 
expectation values of quadratic Wilson loops of length r = ka, where a(f5) is 
the lattice spacing and k is an integer. It is convenient for the extrapolation 
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Figure 1: Creutz ratios: fit to data of the full theory (solid line), m- 
projection (circles) and MAG projection (crosses). 

to the continuum limit to introduce the scale A via 

A2 = 012 ^) exp {-TT (fl - 23) }' < 10) 

which is a renormalization group invariant quantity when one-loop scaling 
applies in the asymptotic /5-region. The normalization of A 2 is chosen for 
reproducing the full SU(2) string tension, i.e. A 2 — > a for (3 ^> 2.1. 

Figure [I] shows our numerical data for Xkk m units of A 2 as function of r A. 
It turns out that the data of the simulation employing the full SU(2) Wilson 
action is best fitted by (solid line) 

XkkA 2 = 7i + 12 1 r 2 • (11) 

This ansatz for Xkk is expected by relating Creutz ratios to the derivative of 
the static quark potential. The second term of the latter equation refers to 
the Coulomb interaction while the first term is present due to a non- vanishing 
string tension. Figure [I] also shows the data points for Xkk calculated from 
links Ufj,{x) which are obtained from m-projection (see (0)). These data 
are contrasted with the result for \kk calculated with MAG projected links. 
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Table 1: The quality of projection, i.e., uj (0), as function of /3. 

(3 2.1 2.2 2.3 2.4 2.5 

[m-proj.] 0.932(6) 0.934(4) 0.936(3) 0.938(3) 0.940(0) 
uj [MAG proj.] 0.67(0) 0.68(4) 0.70(2) 0.72(1) 0.73(8) 



In any case, the asymptotic, i.e. rA ^> 1, value for the string tension is 
reproduced within the numerical accuracy. A striking feature of figure [I] is 
that the data from m-projected links agree with the full result for the sizes 
r explored in figure |1[ 

For a quantitative study of the error induced by projection, we calculated uj 
(P) for the case of the m-projection and the MAG projection, respectively (see 
table [TJ). One observes that uj is much bigger for m-projection rather than 
for the case of MAG projection. In particular, the components of the (gauge 
fixed) links U Q which are dropped by projection, i.e. U Q —>■ U, are roughly 
of 5% in size while the generic error due to projection in the case of MAG is 
of order 30%. One therefore expects that not only the string tension but also 
other physical observables which are calculated with m-projected links are 
well described. On one hand this feature is highly desired for constructing 
effective theories which cover a wide span of low energy properties of Yang- 
Mills theory. On the other hand, the introduction of additional degrees of 
freedom (compared with the case of MAG) obscures those ones which are 
responsible for confinement. 



3 Color alignment in m-gauge 
3.1 M-vector correlations 

Let us assume that we are investigating a physical observable which possesses 
a correlation length £ by comparing the full with the m-projected theory. If 
the color vector m(x) is uniquely oriented in a space-time domain of size 
I ^> £, one would recover the standard MAG scenario (provided that I is 
bigger than the scale set by the critical temperature, i.e., I > 0.7 fm, for 
detaining the Casimir effect). In particular, the dual superconductor picture 
would be expected operating if the string tension is the quantity of interest. 
For investigating the existence of such domains, and, in case, for relating 
the infra-red physics in the m-gauge to the well studied physics in MAG, a 
thorough study of the space-time correlations of the vectors m(x) is highly 
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Figure 2: Raw data of the probability distribution of the scalar product rj 
between two neighboring color vectors fh (left) and function f(rj) (right). 



desired. 

The space-time dependence of the gauge transformation Q which maximizes 
the functional Sa x @ induces a correlation of the unit vectors m(x) in space- 
time. For revealing these correlations, we numerically calculated the prob- 
ability distribution of finding a particular scalar product rj of two vectors 
m located at neighboring sites of distance a(/3). The raw data of this dis- 
tribution are shown in the left panel of figure |^. A random distribution of 
vectors would correspond to dP/drj = 1. We clearly observe a maximum of 
the probability distribution at 77 = 1. We find a color ferromagnetic corre- 
lation between the color vectors fh(x). The value at the maximum position 
increases for increasing /3, i.e., for a decreasing distance a{(3) between the 
neighboring vectors. This indicates that m(x) which constitutes a lattice 
vector model so far will become a smooth field in the continuum limit. 

For an interpretation of these results in the scaling limit a(/3) — > 0, it is useful 
to parameterize the probability distribution as follows: 

— = 1 + cexpj —\ f{r)) , (12) 

where the function f{rj) satisfies without a loss of generality the constraints 

C f( V )d V = 0, /(1) = 1- (13) 
Jo 

The crucial finding is that the constants L and c as well as the function f(r]) 
are universal, i.e., independent of the renormalization point specified by (3 
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L=0.12fm 
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P 



Figure 3: The scaling of the color ferromagnetic correlation length L. 

(see figure ^| right panel and figure §). Comparing the numerical data for 
a(P)/L with the asymptotic one-loop /9-dependence (dashed line in figure |3]) 

a(J3) « 0.16fm exp|-^-/3| (14) 

(where a string tension er = (440 MeV) 2 was used as reference scale), the 
extrapolation to the continuum limit yields 

L = 0.1(2) fm, c = 0.50(6) . (15) 

Note that for observing color ferromagnetic domains a probability distribu- 
tion dP/drj(P — > oo) is required which diverges at r\ — 1. However, our 
numerical data obtained in the scaling window (3 G [2.1,2.5] agree with a 
finite value for dP/drj at r\ — 1. Further numerical investigations (e.g. of the 
volume dependence of the distribution) are necessary for a definite conclusion 
on this issue. 

In conclusion of this section, we find color ferromagnetic interactions between 
neighboring color vectors in which increase for decreasing distance a(/3) thus 
indicating that the vector field m(x) is smooth in the continuum limit. In the 
scaling limit, we find that these correlations extend over a range of roughly 
0.12 fm. The color ferromagnetic interaction between the vectors is (most 
likely) not strong enough to induce the formation of color ferromagnetic 
domains in space-time. 
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Figure 4: Quality of projection lu and the color ferromagnetic interaction 
strength /i (left panel); fraction of neighboring vectors with a scalar product 
larger than 0.95 and 0.99, respectively (right panel). 

3.2 Interpolating gauges 



Subsection £]2]has demonstrated that the m-gauge is well adapted for projec- 
tion. Unfortunately, the distribution of the auxiliary color vectors m(x) does 
not support an arrangement of these vectors in domains of constant orienta- 
tion therefore impeding an interpretation of the m-gauge as local realization 
of MAG. For taking full advantage of the elaborated studies of physics in 
MAG ||, 10 1, we generalize the m-gauge condition @ for allowing a smooth 



interpolation between the MAG and the m-gauge by virtue of a gauge fixing 
parameter k. The generalized gauge fixing action [] is 

Sfix = \ £ tr{^(x)m(x) (f/^(x)m(x)j (16) 



H, {x} 

+ « E 

H,{x} 



m(x) m(x + fi) 



maximum , (17) 



Note that the additional term QT7T ) also respects the reflection symmetry 
rh(x) — > —m(x). For k — 0, one recovers the m-gauge (0). For k ^> 1, on the 
other hand, there is a large penalty in action S'fix for non-uniformly oriented 
color vectors m. One therefore obtains the MAG for sufficiently large k. For 
quantifying the color ferromagnetic interaction strength, we introduce 

/i := — — (m(x)m(x + fi)) . (18) 



A^link 



fi,{x} 



1 We thank Torsten Tok for helpful discussions on useful extensions of S'fix 
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Figure 5: Spatial hypercube of one Monte-Carlo sample: region of aligned 
color vectors rh(x) for k = 0.6. 

One finds n — 1/2 for a random distribution of in G RP2, and retrieves the 
MAG for (i=l. 

Figure |] shows our numerical results for the "quality of projection", i.e. u 
(|9|), and as function of k for a 12 4 lattice and for (3 = 2.4. As expected, the 
strength parameter 11 gradually increases with rising k while u monotonically 
decreases. The minimal error by projection is obtained in m-gauge (k = 0). 

The strength parameter fi at large values of k indicate that regions of uni- 
formly oriented color vector rh form. Figure |] also shows the fraction of 
vector pairs m(x), rh(x + //) which possess a scalar product larger than 0.95 
(0.99). The data are obtained on 12 4 lattice and for (3 = 2.4. For illustrating 
the regions of aligned color vectors at large values of k, figure || presents the 
spatial orientation of the color vectors rh for one Monte-Carlo sample at a 
given time slice. The sample was obtained for k = 0.6. A reference vector 
was chosen at the center of the spatial hypercube. If the scalar product of 
a vector rh located at the position x with the reference vector exceeds 0.95, 
an elementary cube which is spanned by the four points x, x + fi, \i = 1 . . . 3 
is marked. One observes that a particular region of (approximately) aligned 
vectors rh is multi-connected and extends all over the lattice universe. This 
property of the regions of alignment does not match with its analog in solid 
state physics, i.e. the Weiss domains of ferromagnetism. 
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4 Conclusions 



In the Maximal Abelian gauge (MAG), a uniquely oriented color vector m 
defines the embedding of the residual U(l) into the SU(2) gauge group. Evi- 



dence has been accumulated || [K| that in this case an (Abelian) dual Meiss- 
ner effect confines particles which carry color-electric charge with respect to 
the U(l) subgroup. Since colored states which are, however, neutral from the 
viewpoint of the U(l) gauge group escape the confining forces provided by 
the dual superconductor mechanism, a refinement, i.e., a non-Abelian ver- 
sion, of the dual Meissner effect is highly desired. The concept of "hidden 
monopoles" is one possibility [|TT|]. 

By generalizing the MAG gauge condition, we have here proposed another 
possibility for a non-Abelian version of the dual Meissner effect. The new 
gauge (m-gauge) admits a space-time dependent embedding, characterized 
by the color vector m(x), of the residual U(l) into SU(2) gauge group. The 
space-time dependence of m{x) is self-consistently chosen to achieve the min- 
imal error induced by projection. It turns out that the color vector m{x) does 
not change under "small" gauge transformations of the link variable. Thus, 
the field rh(x) carries gauge invariant information encoded in the link vari- 
ables. Our numerical results show color ferromagnetic correlations of these 
vectors m which extends over a range of ~ 0.1(2) fm. The strength of these 
correlations seems to be too small for causing the formation of color ferro- 
magnetic domains. 

For relating the m-gauge to the MAG, we have introduced a class of gauges 
which smoothly interpolates between the MAG and the m-gauge by virtue 
of a gauge fixing parameter k. For a wide span of k, the vacuum decomposes 
into multi-connected regions which are characterized by uniquely oriented 
vectors, and which extend all over the lattice universe. The internal structure 
of these regions define an intrinsic lenght scale Iq. Each region bears the 
potential of an Abelian Meissner effect which operates with respect to the 
residual U(l) subgroup of SU(2). Colored states which do not feel a confining 
force in one particular region generically carry charge in another sector of 
space time. We speculate that, on performing the Monte-Carlo sampling, all 
colored states are confined on length scales bigger than the intrinsic size Iq 
of the regions of color alignment. Note that the size Iq is controlled by the 
gauge parameter k. The request that the average size is a physical quantity 
defines the "running" of the gauge parameter, i.e., the function k((3). The 
actual size of the regions of color alignment in physical units then defines the 
renormalized value Kr and must be provided by a renormalization condition. 
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In subsumption, for a class of gauge conditions, specified by a gauge param- 
eter k > 0, the vacuum consists of regions of aligned color vectors m(x). The 
m-gauge appears as the limiting case k — 0. In this case, the error induced 
by projection is minimal at the expense of additional degrees of freedom as 
compared with the MAG. This fact renders the identification of the degrees 
of freedom relevant for quark confinement more difficult than in the MAG, 
but makes the m-gauge a convenient starting point for formulating an effec- 
tive theory covering a wide span of low energy properties of SU(2) Yang-Mills 
theory. 
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